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Abstract. Using asymptotic Pade approximant methods, we have calculated the 0(0;^) contri- 
bution to the Standard-Model H — > gg decay rate. This process is of particular interest because 
^ \ of the slow convergence evident from the known terms of its QCD series. The 0(as) term is 

expressed as a 3rd degree polynomial in L = In ^ ^ • We find that the asymptotic Pade pre- 
dictions for the renormalization-group accessible coefficients of L, and are within 1%, 2% 

■ and 7% of their respective values. The inclusion of the ©(a^) contribution renders the full decay 

■ rate virtually scale independent over the entire Q.ZMh < H < Mt range of the renormalization 
scale ^. 

'ph: introduction 

D I The Higgs boson remains the only missing particle within the Standard Model 

'-^ I of particle physics. One of the principal hadronic decay modes for a Weinberg- 

Salam Higgs (of mass 100 — 175 GeV) is the Higgs — > 2 gluon {H gg) process, a 
^ ■ process whose leading (one-loop) contribution is 0{al). The two known subleading 

^ , contributions exhibit very slow convergence: e.g. the decay of an Mh = 100 GeV 

Standard-Model Higgs is characterized by the following perturbative series[l]: 

1 + 0.66 + 0.21. 

The present paper is directed towards an estimate of the next order (i.e., four- 
loop) contribution of the H ^ gg process by using renormalization group (RG) 
and asymptotic Pade-approximant methods similar to those employed in [2] for 
next-order terms in correlation functions, and in [3] for the case of a (non-Standard- 
Model) CP-odd Higgs decay into two gluons. In the analysis of this latter process, 
we obtained Pade-approximant estimates of RG- accessible coefficients of next-order 
perturbation theory that are accurate to within a few percent, providing support 
for the Pade-predicted RG-inaccessible coefficient still needed for an estimate of 
the 0{ag) contribution. In the present work, we derive the RG accessible coef- 
ficients of the Standard-Model H ^ gg rate within an MS framework with six 



active flavours [1]. RG invariance aUows us to calculate all the coefficients Ck of 
except for the constant {k — 0) term. We then demonstrate how Fade 



approximant methods may be utilized to estimate the full 0{a\) contribution to 
the decay rate, and show that such estimates of the RG-accessible coefficients are 
quite close to their true values. Finally, we discuss how incorporation of the es- 
timated four-loop term eliminates virtually all residual scale- {fi-) dependence in 
the overall H ^ gg rate. 



RG-ACCESSIBLE COEFFICIENTS OF THE H ^ gg 

RATE 

The H ^ gg decay rate is given by the following perturbative expression in 
the » Mjj» Ml limit [1,4]: 
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The rate (1) is expressed entirely in terms of a scale-dependent i-quark mass 
and coupling constant characterised by six active flavours via L(iji) = In (^^^r^) 



and x^^\iJ,) — a'f\n)/7r. The constant T = In involves a scale-independent 

ratio of propagator-pole Higgs and t-quark masses [4]. The coefficients q in (2) 
characterize the unknown 4th loop contribution to the decay rate. Three of these 
unknown coefficients [ci, C2, C3] can be extracted via the RG-invariance of the phys- 
ical decay rate: 







' ^a;,L,T], 



= [i-2t™W]5^ + /?W^. 

The P and jm functions for six active flavours are given by 
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^^^\x) = -X -'^^x^.. . (5) 

One can verify that (3) is valid order by order to known terms in (2). The continued 
vahdity to orders x^L"^, x^L and x^ yield values for ci, C2 and C3: 

16643^ 3703_, 1225 1127_ 343 

= 910.3167 - —T + —T\ c. = — - — T. = (6) 

In the following section, we will utilize asymptotic Pade approximation procedure 
to estimate {cq, ci, C2, C3}. The ability to predict known coefficients will serve to 
test the appropriateness of this procedure for estimating the RG-inaccessible coef- 
ficient Cq. 



PADE PREDICTIONS FOR FOUR-LOOP 
COEFFICIENTS 

The series S[x,L,T] in (2) may be expressed as follows: 

S[x,L,T] =x^ [l + Ri[L,T]x + R2[L,T]x^ + R^[L,T]x'^ + 

where 
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Rs[L, T] = co(T) + ci(r)L + C2(r)L2 + c^L^ 



(9) 
(10) 



As discussed in prior work [2, 3, 5] and in another contribution to these Pro- 
ceedings,[] Pade-approximant methods may be utilized in order to predict that 



Rs[L] = 2Rl[L]/{R,[L]R2[L] + i??[L]) 



(11) 



To obtain ci, 02,03, we match the scale dependence of (10) to that of (11) over 
the L > region [ultraviolet scales of yU > Mt(/i)] by optimizing the least squares 
function 



X^[co, Ci, C2, C3] 
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dw, 



(12) 



as in our prior contribution to this Proceedings,^ with w = ^^^2 [L = —ln{w)], 
and with R3 given by (11). For Mh = lOOGeV (T = -1.126), we find that 

X^(co,Ci,C2,C3) 



^See "Reducing the Theoretical Uncertainty in the Extraction of \Vub\ ■ ■ •" in this vohime. 



= 4.064168878 • 10^ + 720c^ + 4coC2 + + 2coCi + 24c^ + 12ciC2 + 12coC3 + 20^ 
+ 24OC2C3 + 48ciC3 - 10205.87co - 17507.05ci - 54106.23c2 - 234536.67c3. (13) 
Optimization requires that 



which yields the following predicted values: 

ct = 2452, = 1774, ct = 377.2, c^' = 20.45, (15) 

The corresponding true values of the RG-accessible coefficients (6) when — 
100 GeV are 

ci = 1789.02, C2 = 385.568, C3 = 21.4375. (16) 

indicative of the accuracy of the estimation procedure leading to (15). In Table 1, 
a tabulated set of predictions of the four loop term Cq is given for Mu between 100 
and 175 GeV. Relative errors of Ci and C2 are seen to remain within 2% over the 
range of Mh given; C3 is seen to be within 7% of its true value for the same range. 

The results for cq, as tabulated in Table 1, can be utihzed to predict cq as a 
degree-3 polynomial in T: 

co(T) = 755.9 - 1029T + 394.3T2 - 26.74T^ (17) 

This prediction, as obtained via Pade-approximant methods, can be coupled with 
C3 and the known T-dependence (6) of Ci(T) and C2(T) to predict the full four-loop 
contribution (10) to S[x{iJ,), L{iJ,),T], the scale-sensitive (i.e., /^-dependent) portion 
(7) of the H ^ gg rate (1). 

Surprisingly, the scale-dependence of S[x{ii), L{iJ,),T] is virtually eliminated 
upon incorporation of the estimated four-loop contribution to the rate. In Table 2, 
this scale-sensitive portion of the H ^ gg rate is evaluated for = 125 GeV for 
selected values of the scale parameter between 38 and 175 GeV. The four-loop 
term R^x^, as defined by (10) to include the Mh — 125 GeV Fade estimate Cq = 
1646 from Table 1, varies from zero to 7% of the leading term in the perturbative 
series within S[x,L,T], depending on the choice for /i. Nevertheless, S[x,L,T] 
displays a relative error of 0.6% over this same range of fi, suggesting that the 
Fade estimate for cq is precisely what is required to remove virtually all residual 
scale- (/X -) dependence from the four-loop-order decay rate. 
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TABLE 1: Predicted (pr) and true (RG) values (as determined by renormalization- group meth- 
ods) for four-loop coefficients of the H gg decay rate, as obtained via the moment approach 
of ref. [4]. Least- squares predictions, as in (15) for Mh — 100 GeF are virtually identical to 
those of the moment approach, which also depend on the asymptotic Pade prediction (11) for 
the four-loop term. 
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1+ 0.7127 + 0.2668 + 0.0676 


0.002475 



TABLE 2: Virtual scale invariance oi H ^ gg four-loop estimated rate when Mh = 125 GeV. 

as{n) is assumed to evolve from a.s{Mz) = 0.119 [6], and mt^n) is assumed to evolve from 
mt{mt) = Mi = 175.6 GeV^. The order-by-order perturbation series within the decay rate, as 
defined by (7), as well as S[x{ijl), L{iJi), T], the scale-sensitive portion of the rate (1), are listed for 
representative values of fi between Mh/3 and M^. Although the predicted four loop term Rsx^ 
varies between zero and seven percent of the lead (one-loop) term in the series, the scale-sensitive 
portion of the rate S[x, L.T] is constant up to | AS/S'] = 0.6%. 
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TABLE 3: Virtual scale invariance oi H ^ gg four-loop estimated rate when Mr = 150 GeV. 
Running of as{n) and mt{fJ,), as well as the perturbative series I+R1X+. . . within S[x{iJ,), L{ij,), T], 
the scale-sensitive portion of the rate (1), are as in Table 2. Although the four-loop term Rax^ 
grows to be 4.6% of the leading term (unity) in the series, the rate itself is seen to exhibit a 
relative error |AS'/S'| < 0.3% over the range Mh/S < fi < Mt. 
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TABLE 4: Incorporation of Mass-corrections 6m{l^^^) to the perturbative series (22) in the 
predicted Standard- Model H ^ gg decay rate for various choices of Mh, as described in the 
text. The "minimal sensitivity" choice of scale (yu^^) represents a (very weak) local maximum 

of S[x{^j,L{fi),T] prior to the incorporation of mass-corrections. The magnitude of the leading 
mass correction S„i is shown to be several times larger than the estimated four-loop correction 
R^x^, although it is generally smaller than the known three-loop correction R2X^. 



This reduction in scale dependence is even more dramatic as Mh increases to 
150 Gey , as is evident from Table 3. Incorporation of the Pade estimate co = 1120 
(Table 1) into R^x^, as defined by (10) [with RG- values (6) for {ci, C2, C3}], leads to 
a four-loop expression for the rate characterised by a relative error |AS'/S'| < 0.3% 
over values of the scale parameter n chosen between 47 and 175 GeV. It must 
be noted that cq's contribution to S[x{ii), L{iJ,),T] is just cqx^{ii), in itself an 
appreciably scale-dependent contribution. The absence of any concomitant scale 
dependence in S, despite its potentially large contribution from an arbitrary choice 
of Co, demonstrates the success of the Pade estimate for cq in controlling potentially 
large residual scale sensitivity in the H ^ gg rate. 

LEADING FERMION MASS CORRECTIONS 

The rate (1) has implicit within its lead term the assumption that << 4M| 
[1]. Physical departures from this assumption, while small compared to the three- 
loop contribution to the H ^ gg rate as calculated in [1], are in fact somewhat 
larger than the four-loop contributions estimated above. These fermion mass cor- 
rections are most easily incorporated by replacing the leading factor of unity in 
the series 1 + R^x + R2X^ + R^x^ with the following departure from the assumed 
ml « M]j « 4Mf hierarchy of masses [4,7]: 

' - ^ UAt + ReAbf + ilmAhf^ 
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We (somewhat arbitrarily) choose the scale /i to be the "minimal- sensitivity" [8] 
maximum (/i*^"^) of S[x{fi), L{fi),T], which is seen to occur at values of fi tabulated 
in Table 4 for various choices of Mh- The values for mb{fJ''^^^) and mt^jj,^^^) utilized 
in (19) and (20) are evolved from refernce values mt{mt) = 175.6 GeV, mb{mb) = 
4.2GeV [6]. The leading mass correction SmifJ^'^^), as obtained from (18-21), 
is tabulated and compared to the corresponding estimate of the four-loop term 
Rsx^ifi^^^) in Table 4. These mass corrections are clearly seen to dominate over 
the four-loop term in the scale-sensitive portion of the H gg rate, which is now 
given by 

^[x(/^^),L(/^^),r] 



= x\fi''''') [1 + 6mir') + i?iX(/i^^^^) + R2x\fl'''') + R3x'ifi'''')\ . (22) 

Predictions for the full Standard Model H gg rate that are inclusive of these 
mass corrections are also listed in Table 4. Although these leading mass corrections 
are two to six times as large as the estimated four-loop correction, they are seen 
to be smaller (in general) than the three-loop corrections R2x'^{fj,^^). Moreover, 
subsequent subleading departures from the ml « Mjj « iM^ mass hierarchy 
are anticipated to be 0{x6m), which should be comparable to or somewhat smaller 
than the four-loop terms hsted in Table 4. Note that |5m(/i*^'^)| is itself only ~ 2% 
of the leading term (i.e., unity) in the perturbative series for the Mh = 125 and 
150 GeV evident from Table 4, thereby providing justification for the 

three-digit accuracy of the H ^ gg predicted rates for these cases. 
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